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Abstract. 

We review the recent proposal by Lattanzi & Valle of the majoron as a suitable warm dark matter candidate. The majoron 
is the Goldstone boson associated to the spontaneous breaking of ungauged lepton number, one of the mechanisms proposed 
to give rise to neutrino masses. The majoron can acquire a mass through quantum gravity effects, and can possibly account 
for the observed dark matter component of the Universe. We present constraints on the majoron lifetime, mass and abundance 
obtained by the analysis of the cosmic microwave background data. We find that, in the case of thermal production, the limits 
for the majoron mass read 0.12 keV < mj < 0.17 keV, and discuss how these limits are modified in the non-thermal 
case. The majoron lifetime T is constrained to be > 250Gyr.. We also apply this results to a given seesaw model for the 
generation of neutrino masses, and find that this constraints the energy scale for the lepton number breaking phase transition 
to be > 10 6 GeV. We thus find that the majoron decaying dark matter (DDM) scenario fits nicely in models where neutrino 
masses arise a la seesaw, and may lead to other possible cosmological implications. 
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Introduction 

Understanding the nature of dark matter (DM) and its 
origin represents one of the longest-standing challenges 
in particle cosmology. We know from cosmological ob- 
servations [1, 2] that only ~ 5% of the Universe en- 
ergy content is accounted for by normal, baryonic mat- 
ter, while the remaining is in the form of dark matter 
(~ 25%) and of a similarly elusive energy component, 
dubbed dark energy (~ 70%). 

Historically, the neutrino was at first seen as the nat- 
ural dark matter candidate, due to its weak interaction 
with ordinary matter [3]. However, it soon became evi- 
dent that the high velocity dispersion of the relativistic 
neutrinos would erase all density perturbations below a 
critical scale of some tens of megaparsecs [4], thus com- 
pletely spoiling the whole process of structure formation. 
This critical scale is called free-streaming length; dark 
matter candidates with a large free streaming length, like 
the neutrino, are classified as Hot Dark Matter (HDM). 
Nowadays, although we know from neutrino oscillation 
experiments that neutrinos do have mass [5], recent cos- 
mological data [6] as well as searches for distortions in 
beta [7] and double beta decay spectra [8] place a strin- 
gent limit on the absolute scale of neutrino mass that 
precludes neutrinos from being viable dark matter can- 
didates [9] and from playing a direct role in structure 
formation. 

Many candidates for the dark matter particle are 
presently under consideration: among the most popu- 
lar, the supersymmetric neutralino, and the Kaluza-Klein 
particles (see [10] and references therein). Most of this 



candidates share the property of being Cold Dark Matter 
(CDM) particles, because their velocity dispersion, and 
consequently their free streaming length, are so small to 
be practically irrelevant for cosmological structure for- 
mation. This avoids the problem of small scale damp- 
ing of HDM models, and in fact CDM models agree 
well with observations down to scales of several Mpc, 
once mildly non-linear effects are taken into account [2]. 
However, it seems that the CDM scenario is unable to 
reproduce the matter distribution at the smallest scales, 
i.e, on Mpc scales and below (see [11] and references 
therein). First, it predicts a number of dwarf galaxies 
much larger than observed. Secondly, numerical simu- 
lations produce DM halos with very high density cores, 
but this cuspiness is not actually observed in real galactic 
cores. It is still unclear if these are shortcomings of the 
model itself, or instead come from our poor understand- 
ing of astrophysical processes important at the scale of 
interest, or even from numerical issues related to the high 
non-linearity of the phenomena under consideration. The 
problem with the CDM scenario is in some sense oppo- 
site with respect to the HDM one: where the latter pre- 
dicts too little power in the small scale fluctuations, the 
former predicts too much of it. In other words: a HDM 
Universe is too smooth with respect to the observed one, 
while a CDM Universe is too clumpy. 

Between the two limiting cases of hot and cold dark 
matter, lies the so-called Warm Dark Matter (WDM). 
Examples of WMD candidates include the sterile neu- 
trino [12] and the light gravitino [13]. The free streaming 
length of WDM particles is in the Mpc range, thus quite 
smaller with respect to the typical HDM value (hence the 



name). This is appealing because it suggest the possibil- 
ity of keeping the successful predictions predictions of 
the CDM scenario at the intermediate and large scales, 
and at the same time alleviating (and hopefully eliminat- 
ing) the small-scale inconsistencies of the model [14]. 
Here we describe our recent proposal of a WDM model 
linking the problem of dark matter with the issue of the 
origin of neutrino masses [15]. 



Majoron dark matter 

If neutrino masses arise from the spontaneous vio- 
lation of ungauged lepton number there must exist a 
pseudoscalar gauge singlet Nambu-Goldstone boson, the 
majoron [16, 17]. This may pick up a mass from non- 
perturbative gravitational effects that explicitly break 
global symmetries [18]. Despite the fact that the ma- 
jorons produced at the corresponding spontaneous L- 
violation phase will decay, mainly to neutrinos, they 
could still provide a sizeable fraction of the dark matter 
in the Universe since their couplings are rather tiny. 

This scenario was first considered in Ref. [19]; how- 
ever, since then there have been important observational 
developments which must be taken into account in or- 
der to assess its viability, most notably the recent cos- 
mological microwave observations from the Wilkinson 
Microwave Anisotropy Probe (WMAP) [1]. 

In the following we then consider the majoron de- 
caying dark matter (DDM) idea in a modified ACDM 
cosmological model in which the dark matter particle 
is identified with the weakly interacting majoron J with 
mass in the keV range. A keV weakly interacting particle 
could provide a sizeable fraction of the critical density 
Pa- = 1.88 x 10~ 29 h 2 g/cm 3 and possibly play an im- 
portant role in structure formation, since the associated 
Jeans mass mj eam ~ m^/nij lies in the relevant range. 

The majoron is however not stable but decays non ra- 
diatively with a small decay rate F. In this DDM sce- 
nario, the anisotropics of the cosmic microwave back- 
ground (CMB) can be used to constrain the lifetime 
T = r and the present abundance Q.j of the majoron; 
here we show that the cosmological constraints on DDM 
majorons not only can be fulfilled but also can easily fit 
into a comprehensive global picture for neutrino mass 
generation with spontaneous violation of lepton number. 
For definiteness here we adopt the very popular possibil- 
ity that neutrino masses arise a la seesaw [20] . 

Majoron abundance. Although majorons could re- 
sult from a phase transition, we first consider them to be 
produced thermally, in equilibrium with photons in the 
early Universe. In this case the majoron abundance nj at 
the present time fo will be, owing to entropy conservation 



and taking into account their finite lifetime: 

nj(t ) = 43/11 nj(t D ) _^ /t 
«y(fo) N D n Y (t D ) 



(1) 



where to is the time of majoron decoupling, and No de- 
notes the number of quantum degrees of freedom at that 
time. The exponential factor accounts for majoron de- 
cay. If T{t D ) > 170GeV, then N D = 427/4 = 106.75 
for the particle content of the standard model, while in 
the context of a supersymmetric extension of the SM, 
there would possibly be, at sufficiently early times, about 
twice that number of degrees of freedom. Finally, in 
thermal equilibrium the majoron to photon ratio / = 
nj{tD)/ny{tD) is equal to 1/2. The present density param- 
eter of majorons is then, using No = 106.75: 



mj 



1.25 keV 



(2) 



Another possibility is that majorons were produced 
already out of equilibrium. In this case there is a range 
of possible models, that we can write generically as 



mj 



1.25keV 



-toh 



(3) 



where the quantity j3 parametrizes our ignorance about 
both the exact production mechanism, and the exact 
value of No- When J3 = 1, we recover the scenario de- 
scribed above, with / = 1/2 and N D = 427 /4. 

Effects of majoron DM on the CMB. Clearly if the 
majoron has to survive as a dark matter particle it must 
be long-lived, T > ?o. However, a more stringent bound 
follows by studying the effect of a finite majoron lifetime 
on the cosmological evolution and in particular on the 
CMB anisotropy spectrum. In the DDM scenario, due 
to particle decays, the dark matter density is decreasing 
faster than in the standard cosmological picture. This 
changes the time t eq of radiation-matter equality. This 
means that, for a fixed £2/, there will be more dark matter 
at early times, and the equality will take place earlier, as 
illustrated in Fig. 1 . The present amount of dark matter is 
Q. DM = 0.25 for both models; T" 1 = 14Gyr in the DDM 
model. Other relevant parameters are (Oy = 2.23 x 10~ 2 
and h = 0.7. The time at which the blue and red lines 
cross is the time of matter-radiation equality; for fixed 
Rdm, it shifts to earlier times as the majoron lifetime 
decreases. 

The time of matter-radiation equality has a direct ef- 
fect on the CMB power spectrum. The gravitational po- 
tentials are decaying during the radiation dominated era; 
this means that photons will receive an energy boost after 
crossing potential wells. This so-called early integrated 
Sachs-Wolfe (EISW) effect ceases when matter comes 
to dominate the Universe, since the potential are con- 
stant during matter domination. The overall effect is to 



increase the power around the first peak of the spectrum 
as the equality moves to later times. 



the majoron dark matter component, while the subscript 
DP denotes the majoron relativistic decay products. 
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FIGURE 1. (Color online) Evolution of abundances in the 
standard (thin lines) and DDM (thick lines) universe scenario: 
blue/short dashed, red/long dashed and black/solid correspond 
to the matter, radiation and A components, respectively. 



On the other hand, since T > to, we expect the majoron 
decays to make the gravitational potentials vary again in 
the late stage of the cosmological evolution. This will 
induce a similar effect to the one described above, only 
affecting larger scales due to the increased horizon size. 
This late integrated Sachs-Wolfe (LISW) effect results 
then in an excess of power at small multipoles. 

Both effects can be used in principle to constrain the 
majoron lifetime and cosmological abundance. In order 
to carry on a quantitative analysis, we have developed 
a modified version of the CAMB code [21], that enables 
to compute the CMB anisotropy spectrum once the ma- 
joron lifetime and abundance are given in addition to the 
standard ACDM model parameters. 

We stress the fact that even if a keV majoron consti- 
tutes a warm dark matter particle, it actually behaves as 
cold dark matter insofar as the calculation of its effect 
on the CMB spectrum is concerned, since CMB mea- 
surements cannot discriminate between cold and warm 
dark matter. The latter behaves differently from a cold 
one on scales smaller than its free-streaming length X/ s . 
For a particle mass in the keV range, we have Xf s ~ 
IMpc which corresponds in the CMB to a multipole 
I ~ few thousands. 

The formalism needed to account for the cosmological 
evolution of an unstable relic and of its light decay prod- 
ucts, has been developed for example in Refs. [22, 23], 
including the modifications in both background quanti- 
ties and perturbation evolution. We report here the nec- 
essary changes to the evolution equations, using the for- 
malism introduced in Ref. [24]. The subscript J denotes 



Background equations 

The equations for the time evolution of the dark matter 
and decay products energy density are: 

pj + 3-pj = -aTpj, (4a) 
a 

p DP +A-p DP = aTpj, (4b) 
a 

where a is the cosmological scale factor, and the dot de- 
notes derivative with respect to conformal time (hence 
the extra a factor on the right hand side). Note that the 
source term for the decay products involves the dark mat- 
ter density, thus effectively coupling the two equations. 

Perturbation equations 

The perturbations in the DM and DP components 
evolve according to the following set of equations (see 
Ref. [24] for the meaning of the symbols): 

Majoron dark matter. 




Decay products. 

8 D p = -l (h + 29 DP ) + -(8j- S DP ) (6a) 
3 r 

e D p = k 2 {^- oop^j - r - 9 DP (6b) 

o DP = (26 DP + h + 6fj) - ^-rkF D p3 - o DP - (6c) 
15 10 r 

k 

Fdpi = xj—t \£Fdp,i-\ - + l)F D p.e+i] 

U+l (6d) 

-F DP /-, £>3. 
r 

where r denotes the ratio of DP to photons energy 
densities, i.e. r = Pop/Py (any fiducial density scaling 
as aT A will actually do the job). The presence of the dark 
matter density perturbation Sj in the right hand side of 
Eq. (6a) again couples the two sets of equations. 

Statistical analysis 

Parametrization. Two distinct mechanisms effective 
at very different times characterize the effect of DDM 



on the CMB. It is therefore convenient to choose a 
parametrization that can take advantage of this fact. In 
particular, the "natural" parametrization (£lj, F) has the 
drawback that both parameters affect the time of matter- 
radiation equality. It is more convenient to define the 
quantity 



where p b is the energy density of baryons, and f ear i y <C 
fo < T. As long as this condition is fulfilled, the value 
of Y does not depend on the particular choice of t ear iy , 
since the ratio pj/pt is asymptotically constant at small 
times. Given that t eq < T we can use the value of Y 
to parametrize the relative abundance of majorons at 
matter-radiation equality. In order to simplify notation 
let us also define F\$ = r/(10 _18 sec _1 ); in this way, 
Tig = 1 corresponds to a lifetime t ~ 30Gyr. 

The advantage of using the parametrization (Y, F) is 
that, when all other parameters are fixed, the time of 
matter-radiation equality is uniquely determined by Y, 
while the magnitude of the LISW effect is largely deter- 
mined by r. 

We show in Fig. 2 how the two physical effects are 
nicely separated in this parametrization. We start from a 
fiducial model with r^g = and Y = 4.7; all other pa- 
rameters are fixed to their WMAP best-fit values. The 
values of Tig and Y are chosen in such a way to give 
Q.jh 2 = 0.10, so that this fiducial model reproduces ex- 
actly the WMAP best-fit. At a larger majoron decay rate 
T 18 = 1.2, i.e., r^ 1 ~ 27Gyr, the LISW effect makes, as 
expected, the power at small multipoles increase, while 
the shape of the spectrum around the first peak does not 
change, since the abundance of matter at early times has 
not changed. Finally, increasing Y by 20% the height of 
the first peak decreases accordingly, while the largest an- 
gular scales (small is) are nearly unaffected. A small de- 
crease in power in this region is actually observed, and 
can be explained by noticing that increasing the matter 
content we delay the onset of the A dominated era, re- 
ducing the A contribution to the LISW effect. 

Another advantage of using the above parametrization 
is that Y is directly related to the majoron mass through: 

F = 0.71xP) (P). (8) 



long as the present age of the Universe, is quite at vari- 
ance with respect to the WMAP data. However one must 
take into account the fact the values of the other cosmo- 
logical parameters can be arranged in such a way as to 
reduce or even cancel the conflict with observation, i. e. 
degeneracies may be present in parameter space. In or- 
der to obtain reliable constraints for the majoron mass 
and lifetime, we perform a statistical analysis allowing 
for the variation of all parameters. This is better accom- 
plished using a Markov chain Monte Carlo approach; we 
used to this purpose the widely known COSMOMC code 
[25]. 

In our modified flat (£2 = 1) ACDM model, all the 
dark matter is composed of majorons. This means that no 
stable cold dark matter is present 1 . The 7-dimensional 
parameter space we explore therefore includes the two 
parameters (Y, F) defined above, in addition to the five 
standard parameters, namely: the baryon density £2/,/i 2 , 
the dimensionless Hubble constant h, the reionization op- 
tical depth X re , the amplitude A s and spectral index n s 
of the primordial density fluctuations. The cosmologi- 
cal constant density £2a depends on the values of the 
other parameters due to the flatness condition. We com- 
pare our results with the CMB anisotropics observed by 
the WMAP experiment. Once the full probability distri- 
bution function for the seven base parameters has been 
obtained in this way, the probability densities for de- 
rived parameters, such as the majoron mass mj, is conse- 
quently calculated. 

We show our result in Fig. 3, where we give the 68% 
and 95% confidence contours in the (mj, F) plane, for 
the case J3 = 1, i.e., thermal majoron production and 
No = 427/4. We note that these parameters are not de- 
generate one with the other, so the respective constraints 
are independent. Similarly, we find no degeneracy be- 
tween r and the five standard parameters. The marginal- 
ized 1 -dimensional limits for F and mj are: 

r< 1.3 x lO^sec" 1 (9) 
0.12keV<m y <0.17keV (10) 

Expressed in terms of the majoron lifetime our result im- 
plies t > 250 Gyr, nearly a factor 20 improvement with 
respect to the naive limit t > ?o — 14 Gyr, illustrating 
the power of CMB observations in constraining particle 
physics scenarios. 



Results and discussion 

We are now ready to compute the constraints that 
CMB observations put on the majoron abundance and 
lifetime. As seen from Fig. 2, even a lifetime twice as 



Let us comment on the possibility that /3 ^ 1 . From eq. 
3, it can be seen that this amounts to the transformation 
m j — > j3my. For example, as we have already pointed out, 



This happens e. g. in models where supersymmetry with broken R 
parity is the origin of neutrino mass [26, 27]. 
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FIGURE 2. (Color online) Effect of DDM parameters on the CMB anisotropy spectrum. The value of the parameters are as 
follows. Red/solid: fiducial model (T l& , Y) = (0, 4.7). Green/dashed: (r 18 , Y) = (1.2, 4.7). Blue/dotted: (Tig, Y) = (1.2, 5.6). See 
text. 
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FIGURE 3. (Color online) Contours of the 68% (green/dark) 
and 95% (yellow/light) confidence regions in the (Tj,mj) 
plane. 



Nd can be as large as 427/2 ~ 200, so that the above 
limit would read 0.24keV < mj < 0.34keV. In general, 
if we allow for the possibility of extra degrees of freedom 
in the early Universe, we always have j3 < 1 and then 
m j > 0.12keV . If instead majorons are produced non- 



thermally, one will in general have j8 > 1 . 

Recent analysis of the Lyman-a forest data [28] sug- 
gest that the mass of the warm dark matter particle should 
be larger than at least 1 keV (maybe even an order of 
magnitude more), the exact result depending on the can- 
didate under consideration, on the data used and on the 
analysis pipeline. Taken at face value, these results, com- 
bined with the limits we obtain from the CMB, seem to 
exclude the majoron as a viable dark matter candidate. 
However, care should be taken in naively applying this 
results to the model presented here. First, the limits on 
the mass of the WDM particle have been obtained in 
the case of a stable particle. The effect of the decay on 
the growth of density fluctuations should be taken into 
account to reliably compare the predicted matter power 
spectrum to the observations. Second, the results of the 
Lyman-a forest analysis actually depend on the phase 
space distribution of the particles at the time of decou- 
pling, and then ultimately on the production mechanism. 
On the contrary, our results, when quoted in terms of the 
quantity fimj, are completely independent on the produc- 
tion mechanism. This is due to the fact that, as we com- 
mented above, it is a good approximation to consider that 
the thermal velocities of majorons are negligible, as long 
as only the CMB is concerned. This means that the par- 
ticle mass never enter directly in the perturbation equa- 
tions; instead it only enters indirectly through the back- 
ground quantity Q.j °c fimj, and this should be regarded 
as the quantity that is really constrained by CMB ob- 
servations. A production mechanism resulting in a "sub- 



thermal" (i.e., j3 < 1) majoron abundance, will result in 
the same dark matter energy density being shared be- 
tween a smaller number of particles, and then in a larger 
particle mass. The same result of a larger mass for fixed 
Clj can be achieved if the number of quantum degrees 
of freedom at decoupling is substantially larger than its 
standard model value of 106.75, for example in theories 
with larger gauge groups and representations. Finally, it 
should be reminded that there is still no consensus on 
whether the Lyman-a data can be considered fully reli- 
able, due the various systematics that are involved in the 
analysis pipeline. 

Particle physics model. We now briefly comment on 
the particle physics model. The simplest possibility is 
that neutrino masses arise a la seesaw [20]. In the basis 
V,V C (where v denote ordinary neutrinos, while v £ are 
the SU(2) <g> U(l) singlet "right-handed" neutrinos) the 
full neutrino mass matrix is given as 

/ F 3 V 3 Y V V 2 \ ,,,s 

and involves, in addition to the singlet, also a Higgs 
triplet contribution [29] whose vacuum expectation value 
obeys a "vev seesaw" relation of the type V3V1 ~ v\. The 
Higgs potential combines spontaneous breaking of lep- 
ton number and of the electroweak symmetry. The prop- 
erties of the seesaw majoron and its couplings follow 
from the symmetry properties of the potential and were 
extensively discussed in [17]. Here we assume, in ad- 
dition, that quantum gravity effects [18] produce non- 
renormalizable Planck-mass suppressed terms which ex- 
plicitly break the global lepton number symmetry and 
provide the majoron mass, which we can not reliably 
compute, but we assume that it lies in the cosmologically 
interesting keV range. 

In all of such models the majoron interacts mainly 
with neutrinos, proportionally to their mass [17], leading 
to 

T (^vv)^4. (12) 
mj m% 

The limits obtained above from the WMAP data can be 
used to roughly constrain the lepton number breaking 
scale as v\ > 3 x (l0 6 GeV) 2 , form v ~ leV. 

Future perspectives. The massive majoron has also 
a sub-leading radiative decay mode, J — > 77, making 
our DDM scenario potentially testable through studies 
of the diffuse photon spectrum in the far ultra violet. 
A more extended investigation of these schemes will be 
presented elsewhere including other cosmological data 
such as the large scale structure data from the Sloan Digi- 
tal Sky Survey (SDSS) [30]. In contrast, we do not expect 
the data from upcoming CMB experiments like Planck 



to substantially improve our bounds on the majoron de- 
cay rate, since they mainly affect the large angular scales 
where the error bars have already reached the limit given 
by cosmic variance. We also note that direct detection 
of a keV majoron is possible in a suitable underground 
experiment [31]. 

Note added in the arXiv version. After we sent this 
paper for publication, another work appeared, studying 
the possibility that the entropy production and bulk vis- 
cosity associated to the late decay of the dark matter par- 
ticle could explain the observed acceleration of the Uni- 
verse [32]. However the decay rate needed to accomplish 
this seems to be too large with respect to the limit found 
in the present analysis of the CMB data. 
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